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We propose a direct method of detection of the nuclear anapole 
moment. It is based on the existence of a linear Stark shift for 
alkali atoms in their ground state perturbed by a quadrupolar 
interaction of uniaxial symmetry around a direction n and a 
magnetic field. This shift is characterized by the T-even pseu- 
doscalar (n • B) (n A E ■ B) / B 2 . It involves on the one hand the 
anisotropy of the hyperfine interaction induced by the quadrupo- 
lar interaction and, on the other, the static electric dipole mo- 
ment arising from electroweak interactions inside the nucleus. 
The case of ground state Cs atoms trapped in a uniaxial (hep) 
phase of solid 4 He is examined. From an explicit evaluation of 
both the hyperfine structure anisotropy and the static dipole de- 
duced from recent empirical data about the Cs nuclear anapole 
moment, we predict the Stark shift. It is three times the exper- 
imental upper bound to be set on the T-odd Stark shift of free 
Cs atoms in order to improve the present limit on the electron 
EDM. 
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Introduction 



It has been well demonstrated that parity violation in atomic transitions can be used 
to test electroweak theory |IJ. In this way, the Standard Model has been confirmed 
convincingly in the domain of low energies. At present, refinements in experiments 
and theory allow more precise measurements to look for a breakdown of the Standard 
Model predictions and hence, new physics [0, 0, £|, |!|. The essential parameter ex- 
tracted from atomic parity violation (PV) measurements is the weak nuclear charge 
Qw- This electroweak parameter appears in the definition of the dominant electron- 
nucleus PV potential induced by a Z exchange: 

Vff{r) = G F /V2 ■ Q w /2 • 75 • P v (r) , (1) 

where the Z couples to the nucleus as a vector particle, just as the photon does 
in the Coulomb interaction. In this Z exchange, Qw plays the same role as the 
electric charge in the Coulomb interaction. 75 is the Dirac matrix which reduces 
to the electron helicity, a ■ p/m e c, in the non-relativistic limit. The distribution 
P\/( r ) normalized to unity represents the weak charge distribution inside the nucleus. 
The physical quantity measured in atomic PV experiments is a transition electric 
dipole moment, Ef v between states with the same parity, like the nSi/2 — > n'Si/2 
transitions. In particular the 65*1/2 — * 7Si/2 transition in cesium has been the subject 
of several experiments, the accuracy of which has been steadily increasing with time 
§0, 



On top of this the PV electron-nucleus interaction involves also a nuclear spin- 
dependent contribution which can provide valuable and original information regarding 
Nuclear Physics. It is generated by an interaction of the current-current type with a 
vector coupling for the electron and an axial coupling for the nucleus. The associated 
PV potential Vffl is given by the following expression: 

Vff = G F /V2 ■ A w / (21) -a- 1- P A {r) , (2) 

where a is the Dirac matrix associated with the electron velocity operator, I the 
nuclear spin and Pa{ t ) a nuclear spin distribution normalized to unity. The weak 
axial moment of the nucleus, Aw, receives several contributions. The most obvious 
one comes from the weak neutral vector boson Z$ with axial coupling to the nucleons. 
However, in the standard electroweak model the coupling constants involved nearly 
cancel accidentally. As first pointed out by Flambaum et al. ||11|| , a sizeable contri- 



bution to Ay/ is induced by the contamination of the atom by the PV interactions 
between the nucleons which take place inside the nucleus. The concept relevant to 
describe this interaction is the nuclear anapole moment In fact the interac- 



tion can be interpreted simply in terms of a chiral contribution to the nuclear spin 
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Figure 1: Simplified representation of the nuclear helimagnetism (figure adapted from C. Bouchiat 
p3[). The normal spin magnetization M. s (f) is assumed to be a constant vector parallel to the 
nuclear spin, distributed uniformely inside a sphere. Under the influence of PV nuclear forces, 
the nuclear magnetization distribution inside the nucleus acquires a chiral parity non-conserving 
component M pnc (r), obtained by rotating Ms(r) through the very small angle (3{r) around r. Three 
chiral magnetization lines in the equatorial plane are shown. The vertical normal magnetization is 
actually larger than A4 pnc {r) by about six orders of magnitude. It can be shown Jl^j that the vertical 
anapole moment is given, within a constant, by the magnetic moment obtained by identifying the 
chiral magnetization lines with lines of electric currents. 



magnetization [|T^, [14]], as illustrated in Fig 1. In other words, one can say that the 
PV nuclear forces inside any stable nucleus are responsible for the nuclear anapole 
moment or equivalently a nuclear helimagnetism. The present paper addresses the 
problem of how to detect directly this unique static nuclear property characteristic of 
parity violation in stable nuclei. 

Up to now there has been only one experimental demonstration of the nuclear 
anapole moment, namely that obtained very recently by the Boulder group [[II]]. In 
their experiment which gives a high precision determination of parity violation in the 
atomic 6Si/ 2 — > 7Si/ 2 Cs transition, this effect appears as a small relative difference, 
actually ~ 5%, between the E\ v transition dipole amplitudes measured on two dif- 
ferent hyperfme lines belonging to that same transition. In this case the dominant 
source of P without T violation comes from the electron-nucleon Z exchange asso- 
ciated with the weak charge Qw of the nucleus. This makes the extraction of the 
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nuclear spin-dependent part a most delicate matter. In view of the importance of this 
result for the determination of the PV pion-nucleon coupling constant, fl (see |0|), 
a totally independent determination is highly desirable. 

It is well known that T reversal invariance forbids the manifestation of V$ in 
an atomic stationary state. However, we shall show in the following sections that in 
such a state T reversal invariance does not forbid the manifestation of hence 
that of the nuclear helimagnetism. For a free atom, the rotation symmetry of the 
Hamiltonian leads to an exact cancellation of the diagonal matrix elements. This 
property still holds true if the rotation symmetry is broken by the application of 
static uniform electric and magnetic fields. However, if the symmetry is broken by 
the application of a static potential of quadrupolar symmetry, for instance by trapping 
the atoms inside a crystal of hexagonal symmetry, then, the stationary atomic states 
are endowed with a permanent electric dipole moment which can give rise to a linear 
Stark shift. This offers a novel possibility of detecting the nuclear helimagnetism 
having a twofold advantage: 

i) In a stationary state it is the sole cause of P without T violation. 

ii) It manifests itself by a modification of the atomic transition frequencies in 
an applied electric field, i.e. a linear Stark shift, providing for the first time an 
opportunity for demonstrating the static character of this unusual nuclear property. 

There exist in the literature other proposals for a direct detection in atoms of the 
nuclear spin-dependent effect, i.e. without any participation from V^: 

i) One is based on the difference between the selection rules of the potentials Vffl 
and V@'. While the former acts as a scalar in the total angular momentum space and 
mixes only states of identical angular momentum (and opposite parity), the latter acts 
like a vector and mixes states of different total angular momentum. Consequently, 
one can find atomic transitions between states of the same parity which are allowed 
for the nuclear spin dependent contribution but remain forbidden for the nuclear spin 
independent one [pj. One such example is the (6p 2 ) 3 Po — ► (6p 2 ) 1 5o lead transition at 
339.4 nm, strictly forbidden for even isotopes, which acquires a non- vanishing matrix 
element E\ v in odd isotopes owing to the PV interaction involving the nuclear spin, 
which mixes the (ftp 2 ) 1 So state to the (6p7s) 3 Pi state of opposite parity |16| . 

ii) A second approach, invoked by several groups in the past and now under seri- 
ous consideration [17]], consists in the detection of an Ef v amplitude via a right-left 



asymmetry appearing in hfs transition probabilities for the ground state of potas- 
sium in the presence of a strong magnetic field (magnetic and hyperfine splittings of 
comparable magnitude). 

iii) There is also the possibility of detecting the energy difference in the NMR 
spectrum of enantiomer molecules JTS[ . 

In view of the extreme difficulty of these other projects, we believe that, over 
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and above its intrisic scientific interest, the linear Stark shift discussed in this paper 
deserves careful consideration. 

The first section of this paper recalls the main angular momentum properties of 
the permanent nuclear spin-dependent PV electric dipole operator arising from the 
nuclear anapole moment. In addition we compute its magnitude for the cesium atom 
using recent empirical data relative to the Cs 65* — ► 7S transition. The next section 
(sec. 2) shows that this dipole can manifest itself via a linear Stark shift only if 
the free atom symmetry is broken. After this we consider the case where the atom 
is perturbed by a crystal field of uniaxial symmetry. Here, the crystal axis n, and 
the applied electric and magnetic fields create a chiral environment permitting the 
existence of a linear Stark shift, the explicit expression for which is given. In the 
section 3, we examine a realistic experimental situation where its observation looks 
reasonably feasible: this deals with Cs atoms trapped inside a 4 He crystal matrix 
of hexagonal symmetry. We have investigated quantitatively how, by breaking the 
atomic symmetry, the matrix induced perturbation manages to generate a linear Stark 
shift. Moreover, we evaluate both the matrix induced anisotropy and the shift. The 
details of the necessary calculation based on a semi-empirical method are given in the 
Appendix. In the final section we suggest another experimental approach in which 
the atoms are no longer submitted to a crystal field, but are instead perturbed by an 
intense nonresonant radiation field. 

1 The permanent nuclear spin-dependent PV electric dipole 

1.1 Symmetry considerations 

The space-time symmetry properties of the atomic electric dipole induced by the 
nuclear spin dependent PV interaction have been presented before in many review 
papers (see for instance |19j). We recall them here for completeness, since they con- 
stitute the starting point of the linear Stark shift calculation developed in the present 
paper. 

First, we wish to stress that the existence of the anapole moment interaction not 
only implies the existence of a transition dipole proportional to the nuclear spin, 
but also that of an electric dipole operator having diagonal matrix elements between 
stationary atomic states. This electric dipole is found to be proportional to the 
operator s A I. Therefore it does not undergo the same transformation under P as 
does an ordinary dipole, since it is a pseudovector instead of a vector. We also note 
that it is even under T-reversal, so that the quantity (s A I) ■ E, associated with a 
linear Stark shift, violates P, but does not violate T invariance. 

It is convenient to define d pv (n',n) as the effective pv electric dipole moment 
operator acting in the tensor product £$ <8> £i of the electronic and nuclear angular 
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momentum spaces, which describes the transition between two S1/2 subspaces corre- 
sponding to given radial quantum numbers n and n' . This effective dipole operator 
includes both contributions from potentials V^' and Vffi. Rotation invariance to- 
gether with the fact that V&> is linear in / implies that d pv (n',n) can be written 
under the following general form: 

d pv (n, n) = —UmE^ pv (n, n) a + i a(n, n) 1 + b(n, n) s A I , (3) 

where the real quantities a(n, n') and b(n, n') parametrize the contribution of the 
nuclear spin-dependent pv potential. Time reversal invariance of V p ty and Vffi implies 
the following relations under the exchange n n': 

ImEf p \(n,n') = —ImE^ pv (n',n), 
a(n,n') = —a(n',n), 

b(n,n') = b(n',n). (4) 

The effective pv static dipole moment V pv = d pv (6,6) relative to the ground state is 
then given by : 

V pv = 6(6,6) s A I = diS A J. (5) 

If we introduce the total angular momentum F = s + I, using simple relations of 
angular momentum algebra, one can derive the useful identity: 

sAl=[F\^s\. (6) 

It then becomes obvious that, in low magnetic fields and without external pertur- 
bation, the dipole operator T> pv has no diagonal matrix elements between atomic 
eigenstates. In fact, as demonstrated in the next section of this paper, a manifesta- 
tion of this dipole requires special conditions for breaking the free-atom rotational 
symmetry. 



1.2 Magnitude of the permanent dipole. 

The magnitude, dj, of the permanent dipole will play a decisive role in the assessment 
of the feasibility of an experiment. We are now going to perform the evaluation of 
dj in the interesting case of cesium. We proceed in two steps : first we compute 
directly 6(6, 7) from experimental data, then we give a theoretical evaluation of the 
ratio 6(6, 6)/6(6, 7). It is convenient to use the notations of ref [0] and fTJ| and to 
rewrite d pv (6, 7) as: 

d pv {6,7) = -tImE p 1 v {QJ)h + r ] j+tr ] 'aAj\ . (7) 



The nuclear spin dependent potential Vffl induces a specific dependence of the pv 
transition dipole on the initial and final hyperfine quantum numbers, F and F'. In 
order to isolate the Vj$ contribution, we are led, following ref [0], to the introduction 
of the reduced amplitudes (Iff 1 '- 

(7S,F'M'\d pv \6S,FM) 
d ^'^ = {F >M>\a\FM) • (8) 

The amplitudes dp F'iVi 7 ]') are tabulated in Table XXII of ref|7|] and reduce to 
—HmEf v (6,7) for vanishing rj and if. The quantity of interest here is the ratio 
Thf — ^43/^3 4 which is given, to second order in 77 and r/, by : 

21+1 , 

r hf ~l — 77. (9) 



Using the empirical value for the ratio r^f given by the last Boulder experiment [ID 
r hf -l = (4.9 ± 0.7) x 1(T 2 , we obtain: 

rf = ~ (r hf - 1) = (-2.1 ± 0.3) x 10- 2 . (10) 
lb 

We deduce 6(6, 7) by a simple identification: 

6(6,7) = ImE%° (6,7) j rf = (1.04 ± 0.15x)10~ 13 |e|a , (11) 

where we have used for ImEi V (Q, 7) the empirical value obtained in ref [ ID |: 

ImEf (6,7) = (-0.837 + 0.003) x 10" n |e|a . 

To compute the ratio 6(6, 6)/6(6, 7), we are going to use an approximate relation, 
derived in ref [IB] , which relates the potential Vffl to V^>: 

V p V(rl = K A ^2j. I -Vff(r). (12) 

Here Ka is a constant very close to unity which depends weakly upon the shape of 
the nuclear distributions P\/( r ) and -Pa(^); j is the single electron angular momentum 
and since, as we shall see, only single particle states with j = 1/2 are involved, we 
can write hereafter 2j = a. 

This relation, valid for high Z atoms like cesium, hinges upon the fact that the 
matrix elements (^3/2 1 \ns 1/2) involving p 3 / 2 states are much smaller- by a factor 
2 x 10 -3 - than those which involve pi/ 2 states, {n'p 1 / 2 \V p ^\nsi/2) ■ This is easily 



verified in the one-particle approximation since the radial wave functions at the sur- 
face of the nucleus are very close to Dirac Coulomb wave functions for an unscreened 
charge Z . It is argued in ref |13| that this property remains true, to the level of few 
%, when V p $(r) is replaced by the non local potential Uj$(r,r > ), which describes the 
core polarization effects within the R.P.A. approximation^]. 

The contributions of V® to the effective dipole operator d pv (n,n') are given as 
the sum of the two operators: 

A« = P(riS 1/2 )V®G(E n ,)dP(nS 1/2 ), 

= P{riS 1/2 )dG(E n )V® P(nS 1/2 ), (13) 

where d is the electric dipole operator, G(E n ) = (E n — H at0 m)~ l the Green function 
operator relative to the atomic hamiltonian; P{nSi/2) and P{n'S\/2) stand for the 
projectors upon the subspaces associated with the configurations nS\/ 2 and n'Syz] 
E n and E n > are the corresponding binding energies. It follows immediatly from the 
Wigner-Eckart theorem that the operators and can be written as: 

A« = ih(n, n') a ; B {1) = ik(n, n') a . (14) 

Using now the relation given in equation (|T2| ) and the commutation of a ■ I with the 
pseudoscalar one gets the following expressions for and B^>: 

A m = iK A ^h(n,n>) (a-j)a, 

B^ = l K A ^k(n,n')a (a-j). (15) 
w y 1 J 

We arrive finally at an expression for d pv (n,n') which can be used to compute the 
ratio 6(6,6)/6(6,7): 



d pv (n,n') = i(a + K A -^-j) (h(n,n') + k(n,n')) + 

Qw 1 

a A - -K A ^- (h(n,n')-k(n,n')) . (16) 

1 LJw 



3 To check the validity of the relation ( |l2| ) we have compared the values for 77 and tf obtained in 
this way with those deduced from a direct computation [^o) of d pv (6, 7). The two results agree to 
better than 10%. 
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Time reversal invariance implies h(n,n) = —k(n,n) so that we can write the sought 
for ratio 6(6, 6)/6(6, 7) as: 



6(6,6) 2ft(6,6) 

6(6,7) h(6, 7) - jfe(6, 7) ' 1 J 

The amplitudes h(6, 6) , /i(6, 7) and k(6, 7) can be computed from the formulas given 
in Eqs. (O) and ($M- We have used the explicit values of the radial matrix elements 
(parity mixing and allowed electric dipole amplitudes) for the intermediate states^] 
6-P1/2 — 9Pi/2 and the energy differences involved, which are tabulated in ref. |p0| 
(Table IV)[]. We obtain in this way: 

6(6, 6)/6(6, 7) = 4.152/1.86 = 2.27 . (18) 



Combining the above result with the value of 6(6, 7) given by equation (|TTD we obtain 
the following estimate for df. 

di ~ 2.36 x 10~ 13 |e|a , (19) 

believed to be about 15 % accurate. 

It is of interest to compare dj with the P-odd T-odd EDM of the Cs atom ob- 
tained from a theoretical evaluation using the latest experimental upper bound for 
the electron EDM pT[ : 

|4| < 7.5 x 10- 19 \e\a . 



Using for the cesium anti-screening factor the theoretical value ||22|| : 120 ±10, one gets 
the following upper bound for the cesium EDM, namely the experimental sensitivity 
to be reached for improving the existing bound on |cf e | : 

\dcsEDM I < 9.0 x 10" 17 |e|a . (20) 

We are going to use Eqs.© and (|i~9l) for calculating the linear Stark shift. It is 
interesting to note here that these equations predict also the magnitude of the pv 
transition dipole involved in an eventual Cs project which would be based on the 
observation of hyperfme transitions in the Cs ground state, analogous to the potassium 



project mentioned in the introduction (see also [17]). Therefore both a project of this 
kind and the linear Stark shift discussed here aim at the determination of the same 
physical parameter, dj, but only the observation of a dc Stark shift would prove its 
static character. 



4 We use here the fact that, as noted by several authors, most of the sum (« 98%) comes from 
the four states 6P 1 / 2 , 7P\ii, &P1/2, 9Pi/2- 

5 Note that a misprint in table IV of ref |^0| has caused an interchange between the contents of 
columns 1 and 2 of its lower half (entitled "7S perturbed"). 
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2 The linear Stark shift induced by Vffl 

2. 1 Need for breaking the rotation symmetry of the atomic Hamiltonian 

The parity conserving spin Hamiltonian in presence of a static magnetic field Bq is: 

H sp in = As ■ I — gs^BS ■ B — 7// • B . (21) 

From section 2, we have seen that, to first order in the electric field, the effect of 
Vffl in presence of an applied electric field can be described by the following Stark 
Hamiltonian: 

Hl t v = d I sAl-E = -d I l -[F 2 , s-E]. (22) 

We have noted that the above identity implies the vanishing of the average value of 
Hpl in the low magnetic field limit. We are going to show that this null result still 
remains valid for arbitrary values and orientations of the magnetic field. 

To do this we consider the transformation properties of both H spin and under 
the symmetry B, defined as the product of T reversal by a rotation of 7r around 
the unit vector u = E A B /\EB \, the rotation R(u,ir). It should be stressed that 
the rotation R(u, 7r) and the symmetry O considered here are quantum mechanical 
transformations acting only on the spin states. The external fields are considered as 
real c— numbers and are not affected. One sees immediately that H sp i n is invariant 
under the symmetry 9 = T R{u, ir), while if** changes sign. We conclude that, in 
order to suppress the linear Stark shift cancellation we have to break the G symmetry. 

This symmetry breaking can be achieved, for instance, by perturbing the atomic 
S1/2 state with a crystal field compatible with uniaxial symmetry along the unit vector 
n. A practical realization looks feasible, since it has been demonstrated that Cs atoms 
can be trapped in a solid matrix of helium having an hexagonal symmetry |23|] (see 
also section 3). In this case the alkali S state is perturbed by the Hamiltonian^: 




where both and p = r/a^ are expressed in atomic units. The perturbed atomic 
state is now a mixture of S and D states, with no component of the orbital angular 
momentum along the n axis. The spin Hamiltonian is modified and an anisotropic 
hyperfine interaction is induced by the D state admixture. The new spin Hamiltonian 
reads: 

H sp in = A ± s- 1 + (A\\ - A ± )(s- n){I- n) - g s p B s- B - lit ■ B . (24) 

6 It has been shown |^3| that in the bubble enclosing the cesium atom there is a small overlap be- 
tween the cesium and the helium orbitals. As a consequence, the axially symmetric crystal potential 
inside the bubble can be well approximated by a regular solution of the Laplace equation. 
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It is easily verified that, if n lies in the (B,u) plane, with non-zero components along 

— * 

both B and it, this perturbed atomic Hamiltonian is no longer invariant under the 
transformation 0. 

Another possible method for breaking the symmetry of H sp i n will be presented in 
section 4. 



2.2 Strong magnetic field limit (y s Bo A±, A\\) 

The anisotropy axis is defined as: 

n = cosipz + sinipx. (25) 

Let us consider the nuclear spin Hamiltonian associated with the restriction of 
H sp in to the electronic eigenstate £(n s ,m s ) perturbed by the quadrupolar potential 
H b (n): 

H (m s ) = A ± m s I z + m 8 (A\\ - A ± ) (sin ip cos if; I x + cos 2 ^4) +-f s B m s - r/iB I z 
= m s ljsB + I z (Aj_sin 2 tp + A\\ cos 2 ip — ^—^-) + I X (A\\ — A±) sin ip cos ip] 



B(m e ) is identical to the Hamiltonian seen by an isolated nucleus coupled to an effective 
magnetic field, B e ^(m s ), having the following components: 

B e J> = —(A\\—A±)smipcosip — - 
Bf = 

777 

B e /f = B - (A ± sin 2 i) + A, cos 2 V) — , (26) 

7/ 

or equivalently : 

B eff = B eff sina ^ B eff = ^ geff = geff 

where 

(Ay — A±) sin ip cos ip m s 

tana — 



-77 B + (A± sin 2 ip + Ay cos 2 ?/>) m s 

In other words, the direction z e ^ of B e ^ can be deduced from the z axis by a 
rotation lZ(y, a) by an angle a around the y axis. Hence, the eigenstates of H^jJ^ are 

— * — * 

\m s rfij >, where rhi now stands for the ^-component of the spin / resulting from / 
through the rotation 1Z(y, —a). 

I z — I • z — I • lZ(y, a)z = cos a I z + sin a I x . 
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We can now compute the linear Stark shift associated with the Hamiltonian given 
by Eq. (|2"2"D , supposing the E field directed along the y axis: 



AE st = (m 8 mi\diE s z I x \m s m ; ) 
= djE m s (mi\I x \rfii) 
= diE m s (mi\I ■ TZ(y,a)x\mi) 

= —djE m s mi sin a . (27) 
If we suppose JiBq <C A\\, A± <C |7s|-Bo and \A\\ — A±\ <C A\\ + Aj_, we obtain: 

A\\ ~ A ± 

tan a ~ -r— r : — - sm w cos w wsina, 

U A w + A ±) 

which yields the simplified expression: 

A\\ — A± 

AE st = —diE m s mj ^ — — sin 2ip . (28) 

In this approximation, AE st can be considered as a modification of the hyperfine 
constant linear in the applied electric field. 

In order to show up the transformation properties of AE st , it is useful to express 
this last result in terms of the two fields, E and B, and the unit vector n which defines 
the anisotropy axis: 

(n-B )(n-E A B ) A» - A± . . 

AE st = -2djm s mA ^ / . (29) 

Bq ^|| ' -™-±_ 

From this expression, it is clearly apparent that the linear shift breaks space reflexion 
symmetry but preserves time reversal invariance. It differs from the P and T violating 
linear Stark shift arising from an electron EDM by the fact that it cancels out when 
the quadrupolar anisotropy of the ground state vanishes. It is also obvious from Eq. 
(53) that, in the strong field limit, the size of the Stark shift depends only on the 
orientation of B relative to E and n and not on the strength of the magnetic field. 
Figure 2 represents two mirror- image configurations of the experiment. 



2.3 Limit of low magnetic fields and small anisotropy 

We now consider the limit \A\\ — A±\ <C 7 s I?o <S A±, An. 

The linear Stark shift can be computed by using second order perturbation the- 
ory. H spin is perturbed by both H^ v and H b (n), the latter being responsible for the 

12 



7^ 

n 



Figure 2: Two mirror- image and T-reversal symmetric experimental configurations corresponding 
to opposite values of the pseudoscalar (n • Bo)(n ■ E A Bq)/Bq. 



anisotropy contribution to H spin , i.e. (A\\ — A±)(s- n)(I ■ n). The fields B and E are 
still taken parallel to z and y respectively. We find: 

AE st (F, M) = 2(A|| - A^E cos t^sva^x 

(F M\s z I x + sJ z \F> M')(F' M'\s z I x - s x I z \F M) 



Efm — Ep'M' 



Since the operator s A I is identical to the commutator [F 2 , — |s], we see that only 
the hyperfine states F' ^ F with M' — M ± 1 contribute to the sum. Therefore, 
in the energy denominator we can neglect the Zeeman contribution which is small 
compared to the hyperfine splitting and, in the sum, we can factorize out the energy 
denominator 2(F — I)A\\(I + |). Since F' = F does not contribute, the resulting sum 
can be performed using a closure relation: 

AE st (F,M) = -P ~ A A ± ] d T lE 1 , sin2^(F M\(s z I x + s x I z )(s z I x -s x I z )\F M) . (30) 

Using standard properties of spin 1/2 matrices, we can transform the diagonal matrix 
element above into j{F M\(I% — I Z )\F M). Once taken into account the axial sym- 
metry of the unperturbed atomic state, it still simplifies to |(F M\(I — 3I 2 )\F M). 
We arrive at the final expression: 

AE st (F, M) = k(F, M) A ^ A± djE sin 2^ , (31) 

A\\ + A±_ 

where 

k(F,M) = 2(F - I)(F M|i(f 2 - 3/ 2 2 )|F M)/(2I + 1) . (32) 

The Stark shift coefficients k(F, M) for ^fCs (1=7/2) are listed in Table 1. 

We note that AE st depends on M 2 so the linear Stark shifts of the Zeeman 
splittings E(F,M) — E(F,M — 1) have opposite signs for M > and M < (see 
figure 3): 

E(4, M) - E(4, M — 1) — hu s M/(2I + 1) + AE st (4, \M\) - AE st (A, \M - 1|) . 



13 




M=-3 ' M=-3 

-2 -2 

-1 -1 

_ 



1 1 



2 2 



3 ' 3 



Figure 3: Energies of the different F, M states of the ground configuration of 1 Cs showing the 
linear Stark shift largely magnified. The two figures correspond to situations realizing opposite signs 
of the pseudoscalar V defined in the text (Left: V > 0; right: V < 0). Top: F=4; bottom: F=3. 
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|M| 





1 


2 


3 


4 


k(4, |M|) 


15/16 


51/64 


3/8 


-21/64 


-21/16 


k(3, |M|) 


-15/16 


- 45/64 





75/64 





Table 1: Linear Stark shift coefficients k(4, M) and k(3, M) of the different F, M substates of the 
natural cesium ground state 



As expected, once again the pseudoscalar 



A\\ -A ± (n ■ B )(n ■ E A B 
A» + A ± Bj 



V = ? ; ^JV^l^I f (33) 



plays an essential role. If V > 0, there is a contraction of the Zeeman splittings 
belonging to the F=4 hyperfine state for positive values of M and a dilatation for 
negative ones, as shown by Figure 3. The situation is reversed when the sign of V 
is changed. In the F=3 hyperfine state, splitting contraction also occurs for M > 
with V > and for M < with V < 0. This behavior could help to discriminate the 
linear Stark shift induced by the nuclear helimagnetism from spurious effects. The 
largest shift between two contiguous sublevels is expected to occur for the couple of 
states F = 3, M = |3| -> F = 3,M = |2|. 

From Table 1 and Eqs. ([19|) and (|31|) we predict: 

75 An — A\ 

AE st (3, 3) - AE st {3, 2) = — -J - ^ sin 2^ d T ■ E , (34) 

with gd 7 ~ 2.76 x 10" 13 |e|ao . 

As in the strong field limit, we note that the size of AE st depends only on the 
direction of B. 



2.4 Analogy between this shift and the PV energy shift searched for in enantiomer 
molecules 

We would like to stress that from the point of view of symmetry considerations there 
exists a close analogy between the linear Stark shift induced by the anapole moment 
and the energy shift which is searched for in enantiomer molecules ||24|| . Indeed in 
the present configuration the three vectors E,B and n which are non-coplanar are 
sufficient to place the atom in a chiral environment similar to that experienced by an 
atomic nucleus inside a chiral molecule. Between two mirror-image environments an 
energy difference is predicted exactly like between two mirror-image molecules. 
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3 Experimental considerations and order of magnitude estimate 

We now consider an experimental situation which looks like a possible candidate for 
the observation of the linear Stark shift discussed in the previous sections. It has been 
demonstrated experimentally p5f that cesium atoms can be trapped in solid matrices 
of 4 He. At low pressures, solid helium cristallizes in an isotropic body-centered cubic 
(bcc) phase, but also in a uniaxial hexagonal close packed (hep) phase. Optically 
detected magnetic resonance has proved to be a sensitive tool for investigating the 
symmetry of the trapping sites. The group of A. Weis has reported the observation 
in the hexagonal phase of zero-field magnetic resonance spectra and magnetic dipole- 
forbidden transitions which they interpret in terms of a quadrupolar distorsion of the 



atomic bubbles ||26|| . Particularly relevant here is their observation of the matrix- 
induced lifting of the Zeeman degeneracies in zero field. This is attributed to the 
combined effect of two interactions, the quadrupolar interaction of the form Hf,(n) = 
Xb Up - n) 2 — | p 2 ^ between the cesium atom and the He matrix on the one hand, 
and the hyperfine interaction in the Cs atom on the other. 

Provided that F 2 is still a good quantum number, it is easily shown from general 
symmetry considerations, that the anisotropy of the hyperfine interaction induced by 
the (hep) crystal potential can be represented, within a given hyperfine multiplet, by 
the effective perturbation: 

H en = C eff {F).{{F.n) 2 - l -F 2 ). 

The constants C e ff(F) can be easily related to the anisotropic hyperfine constants 
appearing in the spin hamiltonian H spin introduced in Eq.(|23|) of the previous section: 



C eff (F = I±l/2) = ±(A ll -A ± )/8. 

In a uniaxial crystal, when the atoms are optically polarized along the crystal axis in 
the absence of external magnetic fields, the lifting of the degeneracy between Zeeman 
sublevels induced by H e jf should make it possible to drive magnetic resonance tran- 
sitions between these levels. One would expect to deduce the hyperfine anisotropy 
from the observed spectra. At first sight, the zero-field magnetic resonance spectra 
observed by Weis et al. would seem to match this prediction. However, their ex- 
periment has been performed in a polycristalline (hep) sample. The effects observed 
in this situation result from averaging over the distribution of the microcrystal axes. 
For each microcrystal, there exists a quantization axis, z, which diagonalizes the hy- 
perfine level density matrix. Immediately a question arises as to the direction of the 
quantization axis z with respect to the microcrystal symmetry axis n. If the popu- 
lation differences resulted, say, from the Boltzmann factor, then z would be along n, 
since in the zero magnetic field limit there is no other preferred direction. In such 
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a situation, there would be no difference between the spectra for a polycrystal and 
a monocrystal. But in the experimental situation considered here, the population 
differences are induced by an optical pumping mechanism which provides a second 
preferred direction: the direction of the photon angular momentum along k. The 
microcrystal density matrix is then expected to keep some memory of the direction 
of k. So, two directions n and k compete in the determination of the quantization 
axis z. To proceed further, we consider the extreme case where z is taken along k, 
together with an assumed isotropic distribution of microcrystal axes. It is then easily 
seen that the lines associated with the hyperfine anisotropy H e ff collapse into a single 
asymmetric line when the average is performed over the polycrystal. Clearly, one at 
least of the two preceeding assumptions is too drastic, most likely the isotropy of the 
n distribution. It is indeed likely that the optical pumping process is more efficient 
for microcrystals having a preferred orientation with respect to the photon angular 
momentum. Such a selection mechanism would then lead to an effective anisotropic 
distribution of n, and in this way a spectrum of separated lines can be recovered. 
From the above qualitative considerations, it clearly follows that the final interpreta- 
tion of the the zero-field resonances requires a detailed analysis of the optical pumping 
process for Cs atoms trapped inside deformed bubbles of arbitrary orientation. The 
corresponding theoretical investigation is currently underway in A. Weis's group. 

Meanwhile, to plan any experiment, we still need to know about the physical origin 
and the magnitude of the ratio , which governs the magnitude of the electroweak 

linear Stark shift. We are going to present now the result of an investigation which 
has led us both to a physical understanding and a reasonably accurate estimate of 
the sought after parameter. We have chosen to devote an appendix to a detailed 
description of our semi-empirical approach, which consists in relating the hyperfine 
anisotropy to another measured physical quantity. Here we shall give a brief summary 
of our procedure and present the final result. 

We start from the remark that there really does exist a mechanism able to generate 
an hyperfine anisotropy to first order in the "bubble" Hamiltonian H^n). The nD 3 / 2 
state is indeed mixed to the 6S1/2 state under the effect of Hb(n), and we note then 
that the hyperfine interaction has non-zero off-diagonal matrix elements between 51/2 
and D3/2 states. In fact, it has been shown previously |27|] that the (nSi/2\Hhf\n'D 3 / 2 ) 
matrix elements are not easy to calculate, because they are dominated by the con- 
tribution coming from many-body effects, due to the existence of an approximate 
selection rule which suppresses the single particle matrix element. However, as we 
show in the appendix, the variation of the matrix elements {n 1 ' S\/2\Hhf\n" 'D3/2) with 
respect to the binding energies £ n 's 1/2 an d £ n "D 3/2 , -expressed in Rydberg- can be 
reasonably well predicted in the limit |£ n 's 1/2 l > l^n" z? 3/2 1 <S 1- In this way, we are 
left with a single parameter which can be deduced from the empirical knowledge of 
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another physical quantity involving the same matrix elements. We have in mind the 
quadrupolar amplitude E 2 induced by the hyperfine interaction which is present in 



the cesium 65 — > 75 transition in the absence of a static electric field |pq| . In order to 
show the relation between the quantities A\\ — A± and E\ , we express them explicitly 
in terms of the matrix elements M(n', n) given by: 

£n'S 1/2 - £n"D 3/2 

The basic formula used in our numerical evaluation of An — Aj_ can be cast in a very 
compact form: 

A A 4 A, 2M(6,6) 

" ~ ~~AS M(7,6) + M(6,7) fl3(7 ' 6) RV (36) 

where AS is the energy of the 65 — > 75 transition and 03 oc / ' \ib is the empirical 
quadrupolar amplitude (see Eq. (A.0) for a precise definition). A second empirical 
input is used to determine the coupling constant \: this is the S—D mixing coefficient 
which is obtained from the hyperfine frequency shifts measured by Weis et al. for Cs 
atoms trapped either in the (bcc) or the (hep) phases p6[ in the low magnetic field 
limit. The ratio involving the matrix elements M(n', n) is evaluated in the appendix, 
using the approximation scheme sketched above. Its absolute value is found to lie 
close to unity. Let us quote now the final result given by our semi-empirical methodF] 
described in the appendix: | A +A± \ = 1-07 x 10 . The uncertainty is believed not to 
exceed 20%. 

For observing the electroweak linear Stark shift discussed in the present paper, it 
is important to work with a uniaxial hexagonal crystal. Indeed, in a polycristalline 
phase, where the individual crystals are oriented totally at random, the average value 
of the pseudoscalar V taken over the isotropic distribution of h is expected to be 
suppressed and thus is the Stark shift computed in the previous section. Although 
trapping of cesium atoms has not yet been achieved in a mono crystalline hexagonal 
phase, the prospect does not look unfeasible and a determination of the magni- 
tude of the hyperfine anisotropy appears to be the first step to be achieved. Using 
- 1.07 x 10 and Eq.(p4[), we find that the effective P-odd T-even electric 



I A||+A_l 

dipole moment of the trapped cesium atoms associated with the nuclear anapole mo- 
ment reaches 2.96 x 10~ 16 |e|ao- For comparison, it is interesting to note that this is 
about three times as large as the Cs EDM limit (Eq.^) to be measured on unper- 
turbed Cs atoms for improving our present knowledge about a possible P-odd T-odd 
EDM of the electron. 



7 This method can be seen as a generalization of that used in sec. 1.2 to evaluate the static dipole 
starting from the empirical knowledge of the transition dipole. 
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4 Breaking the free atom symmetry by application of a nonresonant 
radiation field 

In this last section we want just to mention another possibility for breaking the atomic 
Hamiltonian rotation symmetry by other means than static uniform electric and mag- 
netic fields. We have in mind the application of a strong nonresonant radiation field 
which generates an anisotropic electron gyromagnetic ratio. In the presence of an 
external magnetic field B it has been shown [[HJ that the effect of the nonresonant 
radiation field can be described by the introduction of an effective magnetic field: 

B' = (g ± B + (<7|| - g ± )n ■ B n) /^{g\ + gl) , 

where h defines the direction of polarization of the radiation field, g± = gp and g\\ = 
g^Jo^i/cu), Jo is the zero-order Bessel function, and LO\ is the Rabi angular frequency 
associated with the radiation field. The above formula suggests the existence of a 
uniaxial symmetry, but it is valid only within an atomic hyperfine multiplet. It is 
clear that the "dressing" by a nonresonant radiation field offers new possibilities for 
placing the atoms in a quadrupolar environment. However, it is important to bear in 
mind that at least two stringent requirements must be satisfied if one wants to detect 
an electroweak Stark shift in the ground state. First, the uniaxial perturbation has 
to mix the two hyperfine substates, otherwise the matrix element of if** cancels. 
Second, it is imperative to avoid a broadening of the transition lines for allowing 
precise frequency measurements. We are currently investigating how to achieve the 
proper conditions in a realistic way. 

Conclusion 

This paper investigates a way to get around the well known no-go theorem: 

no linear Stark shift can be observed in a stationary atomic state unless T reversal 

invariance is broken. 

The perturbation of an atom by the nuclear spin-dependent parity-odd potential 
generated by the nuclear anople moment leads to a static electric dipole moment 
di sAl, which clearly is T-even. However, if one considers an atom placed in arbitrarily 
oriented electric and magnetic uniform static fields B and E , the quantum average 
Eq ■ (s A /) is found to vanish. This can be understood by noting that s A I ■ E is 
odd under the quantum symmetry transformation G defined as the product of the 
time reflexion T by a space rotation of tc about an axis normal to a plane parallel to 
the fields B ,E , while the atomic hamiltonian stays even. Our strategy to obtain a 
linear Stark shift is to break the symmetry while keeping T invariance. 

As a possible practical realization of such a situation, we have studied the case of 
ground state Cs atoms trapped in a uniaxial (hep) phase of solid 4 He, which has been 
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recently the subject of detailed spectroscopic studies P6fl . The required breaking of 
space rotation is provided by the uniaxal crystal field. As a result of the deformation 
of the atomic spatial wave function the hyperfine interaction acquires an anisotropic 
part, which plays an essential role in the determination of the size of the linear Stark 
shift. We have performed a numerical estimation of the hyperfine anisotropy, believed 
to be accurate to the 20% level, using a semi-empirical method. We use as an input 
the recent experimental measurement of the E 2 amplitude of the 6Si/ 2 — > 75*1/2 
transition induced in cesium by the hyperfine interaction. We arrive in this way at a 
numerical evaluation of the linear Stark shift induced by the nuclear anapole moment: 
the expected effect is found to be about three times the experimental upper limit to 
be set on the T-odd Stark shift of free Cs atoms for improving the present limit on 
the electron EDM. 

Besides the obvious remark that the T-even Stark shift studied here could be 
a possible source of systematic uncertainty in EDM experiments designed to reach 
unprecedented sensitivity [0, ^T], |32], |33| , we believe that there are strong physical 



motivations for measuring the Stark shift itself. First, it would lead to a direct 
measurement of the nuclear anapole moment in absence of any contribution coming 
from the dominant PV potential due to the weak nuclear charge. It would also provide 
an evidence for a truly static manifestation of the electroweak interaction, something 
which is still lacking. Second, this experiment would rely on the measurement of 
frequency shifts rather than transition amplitudes. While transition probabilities are 
difficult to measure very accurately, high precision measurements of frequency shifts 
have already been achieved. 
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APPENDIX: Semi-empirical calculation of the hyperfine anisotropy of Cs 
atoms trapped inside a 4 He hexagonal matrix 

In this appendix we present our evaluation of the hyperfine structure anisotropy ^— 
resulting from the matrix induced bubble deformation of quadrupolar symmetry, a 
quantity frequently referred to in this paper. 



1. Two processes induced by hyperfine mixing 

Our approach is based on the fact that hyperfine mixing plays quite similar roles 
in two different processes. The first process concerns the Cs 65* — > 7S quadrupolar 
transition amplitude in zero electric field while the second process deals with the 
parameter . 

We start by rewriting the standard mixed Mi — E 2 transition operator in atomic 
units: 

- M 1 

T Ml+E2 = (e A k) hA£(p- e)(p- k) , (1) 

Pb * 

where p is the electron coordinate in Bohr radius unit and AS is the transition 
energy expressed in Rydberg unitf^ We are first going to study the perturbation 
effect on Tm 1 +e 2 caused by the hyperfine interaction Hhf- This phenomenon has 
been observed experimentally in the forbidden 65i/2 — > 7Si/2 transition. It provides 
a useful calibration amplitude in cesium parity violation experiments. To analyse the 
experimental results, it was found convenient, to introduce the effective transition 
operator Thf acting upon the tensor products of the electron spin and nuclear spin 
states: 

T hf = ia 2 (n',n) (s A I) ■ (e A k) +ia 3 {n',n) ((s • k)(I ■ e) + (s • e){T ■ k)) . (2) 

The second physical process to be analysed in this section is not at first sight 
closely connected but happens to be described by the same formalism. This will allow 
us to establish a very useful connection between measurements coming from rather 
different experimental situations. Recently optical pumping has been observed with 



cesium atoms trapped inside an hexagonal matrix of solid helium p6[ . Among the 
new effects to be expected, we have seen earlier in this paper that the existence of an 
anisotropic hyperfine structure opens the possiblity of observing a linear Stark shift 
induced by the nuclear anapole moment, an effect which cannot exist for an atom in 
a spherically symmetric environment. It is known that in the bubble enclosing the 



8 The phase difference, 7r/2, between the two amplitudes expresses the fact that the magnetic 
moment M. and the quadrupole operator behave differently under time reflexion: the first is odd, 
while the second is even. 
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cesium atom there is a small overlap between the cesium and the helium orbitals [23 . 
As a consequence, the axially symmetric crystal potential inside the bubble can be 
well approximated by a regular solution of the Laplace equation: 




The perturbation of the hyperfine interaction by the bubble quadrupole potential 
H b {n) induces an anisotropic hyperfine structure for cesium nS\/2 states. This is 
described by the effective Hamiltonian: 

Hft = ( A W - A ^ (V- ™)( f • *0 - \*' f ) ■ ( 4 ) 

We present now the basic formulas which allow the computation of the parameters 
relevant for the two physical problems in hand. They will be given in such a way as 
to exhibit their close anology. We have chosen to use the Dirac equation formalism. 
Besides the fact that formulas are more compact, it is well known that relativistic 
corrections play an important role in cesium hyperfine structure computation. Ne- 
glecting the contribution of the quadrupole nuclear moment of the Cs nucleus^, the 
hyperfine hamiltonian is written as: 

H hf = I-A, (5) 



A = C hf ^ + 8A^{pJ') + ... . (6) 

The first term gives the hyperfine interaction of the valence electron treated as a Dirac 
particle; the second represents the non-local modification of the hyperfine interaction 
induced by the excitation of core electron-hole pairs to lowest order and the dots 
stand for higher order contributions^]. It has been shown in reference that the off- 



diagonal matrix element (n" D 3 / 2 \^-\nSi/ 2 ) is strongly suppressed by an approximate 

selection rule which does not apply to the many-body non local operator SA^ (p, p'). 
An evaluation of the latter contribution led to a semi-quantative agreement with 
the experimental measurements of the ratio 03(7, 6)/a 2 (7, 6), while the single particle 
result is too small by about two orders of magnitude. 



9 As shown in rcf j28|, the quadrupole contribution for 133 Cs plays a negligible role in the effects 
discussed in this appendix. 

10 An explicit construction of 5A>^(p, p') together with a resummation of an infinite set of higher 
order terms, within the many body field theory formalism, is given in reference p4| . See also ]35| 
for more advanced analysis. 
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To obtain an estimate of the ratio (A\\ — A±) / a${7 , 6) it is convenient to introduce 
the cartesian tensor operator T ili2i3 (E). This object appears naturally in the lowest- 
order pertubation expressions for the quantities of interest: 

T hi2i 3 ( E ) = Ai G£ /2 (E) (pi 2 p i3 - - 5 i2)i3 p 2 ) , (7) 

where the indices i\ , i 2 , 13 take any value between 1 and 3. The scalar operator 
Gy 2 (E) is the atomic Green function operator restricted to the subspace of D 3 / 2 
configurations (total atomic angular momentum J=3/2 and positive parity). We now 
proceed to isolate in T ili2i3 (E) the part transforming as a vector; this is the only part 
to survive after the operator is sandwiched between the projectors P(n'S\/2) and 
P(nS\/ 2 )- This operation is achieved by a decomposition of T ili2i3 (E) into a traceless 
tensor T ili2i3 (E) and a remainder p6| : 



3 

-^1*2*3 (-^0 = -^1*2*3 (-^0 Yq (^1i*2 T a ai 3 \E) + ^Il,l3 ^afCM2 {Ejj 



2 

— 10 ^ i2,is r ^ aaii ^' 



where we have used the fact that T aai = T aia and T iaa = 0. It is a simple matter 
to verify from the above equation that we have indeed T aai3 = T ai2a = T iiaa = 0. 
The fully symmetric part of the traceless tensor T^ li2iz {E) is easily identified with 
an octupole spherical tensor having seven independant components. By a simple 
counting argument, the left over term is seen to have five components; it is to be 
identified with the quadrupole tensor which appears in the full decomposition of 
T ili2i3 (E) into irreducible representations of the rotation group 0(3). Let us have a 
look at the vector operator, V, the components of which appear in the right hand 
side of Eq.(A.|): 

V = (AGt /2 (E) -p)p- l -AG+ /2 {E)p 2 . 

The second term in the above expression does not contribute when it acts upon an 
nS\/2 state so, we are led for our purpose to introduce the vector operator T(n',n) 

f(n', n) = ^P(n'S 1/2 ) (AG+^E^ ■ p) p + (h.c, E f -> ^(^1/2) 

= 7(n',n)s, (9) 

where Ef and Ei are respectively the binding energies of the n' Si/ 2 and nSi/ 2 atomic 
states. The second line of the above equation follows directly from the Wigner-Eckart 
theorem applied to a vector operator. 
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In order to calculate 03(71', n) we have to perform the contraction of I^e^k^ with 
the tensor: 

F ili2 i 3 = P(n'S 1/2 ) (T hi2 i 3 (Ef) + (h.c, E f -> Ei)) P(nS 1/2 ) ■ 
Using Eq.(A|^) and (A|9[) ,F ili2i3 can be cast into the simple form : 



2 

-^1«2«3 liP' 1 n ) ( ^11,12 4~ ^il,l3 $12 2^«2,«3 



The required index contraction with the tensor I^e^k^ is now easily performed and 
one obtains directly a^{n' ', n), up to a prefactor whose value is found by identification 
with Eq.(A.|): 

a 3 (n',n) = -^A£j(n',n) . (10) 

To calculate the hyperfine anisotropy An — A±, we follow the same lines but this time 
the contraction involves the tensor J ia {n i2 n i3 — |5i 2 ,i 3 ), the prefactor is fixed by com- 
parison with Eq.(A|J) and the exchange i 2 H leads to two identical contributions. 
Hence, 

A ]{ -A ± = 2A 6 (^) 7 (n,n) (11) 

e 2 4X b 7(71,71) , 
= - / , r a 3 (n,n). 12 

2ao Ac 7(71,71) 

The expression (A.IT3) looks to us a good starting point for numerical evaluation 
of A\\ — A^: besides the fact that several sources of uncertainties in the evaluation 
of 7(71', n) are eliminated in the ratio ^y" ^ , it lends itself to the use of empirical 
information. One may note, here, a certain similarity with Eq. (7) of sec 1.2 used for 
the evaluation of the permanent dipole, dj. 



2. Numerical evaluation 

We proceed now to a numerical evaluation of A\\ — A± in three steps, starting from 
formula (A.@). 

The numerical value of the 6Si/ 2 — > 7Si/ 2 quadrupole amplitude a 3 (7, 6) is readily 



obtained from measurements 28, of the ratio 



03(7, 6)/a 2 (7, 6) = E 2 /M t l f = (5.3 ± 0.3) x 10~ 2 , 
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combined with a precise theoretical evalution of the magnetic dipole amplitude^ 



a a (7, 6) = = -0.4047 ± 4 x 10" 4 . We obtain finally: 

a 3 (7, 6) = (2.14 ± 0.12) x 10~ 7 . (13) 

The second step is the numerical estimate of the ratio 7(6, 6)/j(7, 6). This is 
more delicate and requires an assumption which has been shown to work in similar 
situations. To begin with, we have addressed the question^ of the origin and size of 
the variations of the off-diagonal matrix elements (nSi/2\Hhf\n" 'D 3 / 2 ) with the radial 
quantum numbers n and n" . It is of interest to remind that a very precise answer to 
this question has been already obtained in the case of cesium single particle matrix 
elements (nLj\H^j\n'Lj) with L = or 1 and with n and n' referring to the radial 
quantum numbers of any pair of valence states. For simplicity, we are going to express 
the answer within a non-relativistic formalism, but it should be borne in mind that 
all of what is said holds true within a relativistic framework. It is convenient to 
introduce the notion of overline matrix elements such as those computed with radial 
wave functions R n ij{p) which have a starting coefficient at the origin equal to unity 
instead of a unit norm0. More explicitly we can write: 

(nLj\H s h p f \n'Lj) 

(nLj\H s h p f \n'Lj) = [ J J f Jl , (14) 

where A ni . = lim p ^ P~ l Rnij{p) is the starting coefficient of the space normalized 
wave function. (In the relativistic case the above condition is replaced by energy 
independent boundary conditions imposed on the Dirac radial wave functions at the 
nuclear radius). It was found in references |27], that the overlined matrix elements 
are independent of the valence orbital radial quantum numbers n and n' to better than 
10~ 4 for S1/2 states and better than 10 -3 for Pi/ 2 states. This result is understood by 
noting that, in the domain of the p values relevant for the evalutation of the matrix 
elements of for S1/2 and Py 2 states, the potential energy is larger than valence 
binding energies by more than three orders of magnitude. This implies that, in this 
domain, the overlined radial wave functions have no dependence upon the binding 
energy or equivalently upon the radial quantum numbers of the valence orbitals. 

11 The theoretical method used to get M± is based upon the factorization rule: {QS\Hhf\7S) = 



hf 



(6S\Hhf\6S)(7S\Hhf\7S). This rule was first established with an accuracy of a few parts in 10 3 in 
ref.^Tj. It has been confirmed by a direct many-body relativistic computation Q of (6S\Hhf\7S), 
accurate to the 1% level. More recently the validity of the rule has been pushed to the level of a 
fraction of 1CT 3 ||. 

12 Arguments similar to those given below and in references ^7), [3?J are developed in |}{| . 
13 The wave function R n ij (p) is known to be an analytic function of the energy. This property is 
the starting point of the quantum defect theory. 
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The above argument has to be reconsidered for the lowest order many body correc- 
tion involving the matrix element of the non local operator: 5A^ (p, p '). The relevant 
domain of p values is now determined by the "radii" of the core outer orbitals involved 
in the computation, which in the case of S\/ 2 and P\/ 2 matrix elements are 5s , 5p , 
while in the case of the off-diagonal matrix element {nSi/ 2 \5A^\n' D 3 / 2 ) only hp is 
relevant. We measure the variation of the overlined matrix elements (nLj\H™j\n' L'j,) 
with the valence state binding energies £ n Lj by the parameters 5lj defined as their 
logarithmic derivative with respect to £ n Lj, (here H^f stands for the many-body mod- 
ification to the hyperfine interaction). From results of references P7| , 57 |, we can infer 

the relative variation of (nLi/ 2 \5A^ \n'Li/ 2 ) for L = 0,1 and we arrive to the values 
8a , = —0.12 and Sp\ = —0.30. The fact that —5p\ is about three times larger 

^1/2 -fl/2 -n./2 & 

than —Sg^ is coming from the fact that P state binding energies have to be com- 
pared with the potential energy minus the centrifugal energy. Let us, now, consider 
the more difficult case of the S-D off-diagonal matrix elements (n' Si/2\SA^\n" -D3/2) . 
The corresponding parameter 5^ is expected to be somewhat larger in absolute 
value, due to the fact that the relevant 5p orbital is less tightly bound than the 5s 
orbital which gives the dominant contribution to the Si/ 2 diagonal matrix element. 
The relative variation versus the D 3 / 2 energy is expected to be on the order of few 
units, since the centrigugal barrier is three times higher than in the case of P states. 
This expectation is borne out by a preliminary estimate which g ives 8j^. ~ — 3. 

We proceed now to a numerical evaluation of the ratio r anis = 7(6, 6)/j(7, 6), 
leaving, for the moment, 5s 1/2 and <5d 3/2 as free parameters. As an intermediate step, 
we compute the quantities M(n',n), written as sums over the intermediate n"D 3 / 2 
states : 

M( n ',n) ~ T ( n ' S ±M H hf\n"D3/2){n"D 3/2 \p 2 \n S 1/2 ) ^ ^ 

n" £n'S 1/2 ~ £n"D 3/2 

The ratio r anis is given in terms of M(n',n) by the following formula : 

r anis = 7(6, 6)/ 7 (7, 6) = ^ g) - ^ f) . (16) 

An explicit numerical computation of r an , s has been performed according to the 
following procedure. First, any binding energy independent factor appearing in 
M(n', n) is dropped since it disappears in the ratio. This is indicated below by the 
symbol oc. The sum J2 n " appearing in M(n', n) is limited to 5 < n" < 8. The set of 
the quadrupole matrix element (n"D 3 / 2 \p 2 \n S\/ 2 ) were obtained by a relativistic ver- 
sion of the Norcross model. In order to test the sensitivity of the result to quadrupole 
amplitudes, we have also used a set calculated by an extension of the Bates-Damgaard 
method. The energy denominators are taken from experiment. The hyperfine matrix 
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elements (n f Si/ 2 \Hhf\n" D 3 / 2 ) , to second order in the binding energies are, given by 
the following formulas: 

(n'Si/ 2 \H hf \n"D 3/2 ) oc A n » D3/2 A n > Sl/2 x 

( X + S S 1/2 £n>S 1/2 + 8D 3/2 £n»D 3/2 ) , (17) 

AnLj cc (S nLj )i. (18) 



In formula (17), we have dropped, according to the above prescription, the zero 



energy limit of the overlined matrix element (p! ' Si/ 2 \Hhf\n" 'D 3 / 2 ) . Equation (A.|18|) 
follows from a result obtained in |[L6|| , where the Fermi-Segre formula was extended 
to arbitrary orbital angular momentum states. For simplicity we have ignored for 
simplicity a factor involving the derivative of the quantum defects, which in the 
present context would introduce few percent corrections. 

We now have all the elements needed to calculate the sought after ratio: 

rams = = -0.8173 - 0.02555 D3/2 + 0.1456<5 5l/2 . (19) 

The negative sign of r anis can be traced back to the fact that the 7Si/ 2 level lies just 
in between 5D 3 / 2 and 6D 3 / 2 levels. If we adopt the rough estimate given above : a 
few tens of % for — 5s 1/2 and a few units for — 5d 3/2 , the first-order energy correction 
remains well below the 10% level, due in part to a cancellation between the two 
correcting terms. To reduce the absolute value of r anis by more than 10% would 
require unrealistic values of — <5d 3/2 so we believe the estimate of r anis = —0.82 ± 0.10 
to be reasonably safe[5 

The final step in our evaluation of the hfs anisotropy is devoted to the empirical 
determination of the coupling constant A& appearing in front of the crystal electronic 
potential. As experimental input we are going to use the hyperfine energy shift which 
is observed for trapped cesium atoms, when one passes from the cubic to the hexagonal 
phase. This shift is attributed to the effect of the anisotropic bubble potential Hb(n). 
We shall ignore, for the moment, the possible contribution of the anisotropic hyperfine 
interaction H%j ts an d assume that the shift is essentially due to the renormalization 
of the 65i/2 component of the atomic wave function by the admixtures a n D 3/2 of the 
nD 3 / 2 states. The corresponding variation of the hyperfine splitting 5W is then given 
by : 

TET = \ a nDj\ 2 = -XlJsD, (20) 

n,J 



14 The validity of the procedure leading to this estimate has been checked, to the 10% level, upon 
a significative subset of the many-body Feynman diagrams contributing to 7(n', n). 
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where we have isolated A 2 by introducing the purely atomic quantity Jsd- Let us write 
down the explicit expression of Jsd, neglecting spin-orbit coupling and assuming that 
n lies along the quantization axis: 



Jsd = E I (651 (cos 2 9 - l/3)p 2 \nD)/(S 6S - S nD )f . (21) 

n 

We limit the sum to n values ranging from 5 to 8. With the same radial quadrupole 
matrix elements as before, we obtain the numerical value: Jsd — 9512. Using the 



empirical number given in ref [26], J—5W/W = 0.035, we arrive at the following 



absolute value of the coupling constant A& (in Ry): 



-SW 

|Atl = Vj^ = a000359 - (22) 

It should be pointed out that if the crystal axis n is not aligned along the quantization 
axis, one obtains values of Jsd smaller than the one quoted above, so the value of 
| Aft | should be considered, strictly speaking, as a lower bound. At last, we have in 
hand all the ingredients needed to perform a numerical evaluation of \A» — A±\ from 
the formula (A.^) since AS = 0.169 is taken directly from experiment: 

\A\\-A ± \ = ^^4^|a 3 (7,6)|Rydberg(MHz)=4.9MHz (23) 
AS |t(7,6)| 

= 1.07 x 10~ 3 x (|A|| + A ± \) . (24) 

As a final topic, we should discuss the effect of the anisotropic hyperfine interaction 
itself on the the empirical splitting SW, since this could modify the value of |A&| 
and so play a role in the assessment of the uncertainty affecting the result given by 
Eq. (A.|4j). Due to this effect, the constant A& is no longer given by Eq. (A.p2j) 
but rather by a second order equation where the linear term is associated with the 
anisotropic hyperfine interaction. It is convenient to introduce the variable x = A&/A° 



with A° = J—SW/(JsdW). The equation giving A& takes then the simple form: 
x 2 — 2b x — 1 =0, where the coefficient b is given by the following formula: 

_ (Ajj - AjjW W A F ( sj z - |£ ■ f) 
A\\+A ± 5W 21 + 1 

The superscript indicates that the hf anisotropy is given, up to a well defined sign, 
by Eq. (Aj24|). The symbol means that one should take the difference between 
the two hyperfine states of the quantum average over which it is applied. To obtain 
an over-estimate of b we have assumed that optical pumping works at its maximum 
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efficiency so that the microwave transition takes place between the hyperfine levels 
(4,4) and (3,3). In this case we obtain b = 0.20 and the two possible solutions for Xb 
are : 

A[ ±} = ±3.610~ 4 (1±0.2). 

The actual experimental situation is expected to lie far from the extreme case consid- 
ered here, so the difference between the two absolute values is certainly smaller than 
the upper limit given by the above calculation. 

In conclusion, including all sources of uncertainties, we consider our evaluation 
of Eq.(A.|24D, | a +a ^ \ = 1-07 x 10~ 3 , as reliable within uncertainty limits of about 
20%. However, if, during hyperfine shift measurements n is not aligned along the 
quantization axis, the central value of A&, and therefore that of , may be pushed 

upwards. 
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